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ABSTRACT 

The study deals with a laathematical model for estimating vegetation cover 
in Tokyo Districts by using LANDSAT MSS data. Two types of model, multi- 
regression analysis and parametric model, were tested for the geographically 
corrected test area with 45 kilometers by 18 kilometers in Tokyo Districts, 
where the ground survey for vegetation cover and classification was carried 
out by the Ministry of Construction, the Government of Japan in 1973. 

The results obtained from the study conclude that vegetation cover can be 
well estimated from LANDSAT MSS data with the accuracy of not less than 80 : 
per cent. r 


i 


INTRODUCTION 

I Vegetation cover can be defined as the ratio of the area with vegetation 

! cover to the mit area. However, recognition of vegetation varies depending 
I on the definition of vegetation. In 1973, the Ministry of Construction carried 
j out the ground survey associated with multi-band camera aerial photographs for 
■ mapping the vegetation and landuse in Tokyo Districts within 50 kilometers 
i radius on the base of 500 meters mesh data. 

j Data for the unit area were digitized with respect to sixteen classes of 

' vegetation and land use, which are utilized as the ground truth data in the 
t study. 
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SELECTION OP TEST AREA 

Test area was selected where Tokyo down town, suburb and hill are properly 
included as shown in Pig. 1. The test area corresponds to two frames of 
national base map, 1:50,000, namely ' HACHTOJI ’ and ' TOKYO SEtNANBU ' .of 
which area makes about 800 square kilometers with 45 km in east*-v 
by 18 km in north-south direction. Since so-^called 500 meters ^ 
by deviding a frame of the base map of 1:50,000 into 40 grids^j 
of longitude and latitude, the test area involves 3200 mesh 
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Pig. 1 


Test Area 



However, 3080 mesli data were utilised because IZO data were not valid for the 
ground survey. 

LAKBSAT MSS digital data were geometrically and geographically corrected 
so as to be registed to the above mensioned base map, of which unit pixel was 
resampled with the size of 100 meters by 100 meters. Therefore, approximetely 
25 pixels ‘ of LAKDSAT data corresponds to a unit 500 meters mesh data* 

GROUND TRUTH DATA 

Vegetation and land use were classified into sixteen classes some of which 
were selected as components of vegetation cover depending on the following three 
cases of definition, as shown in Table 1 , 

(1) Case 1 : Vegetation cover, G^. 

Gi includes ten classes except huilt-'up area, disturbed area, street 
tree, and other. 

( 2 ) Case 2 ; Vegetation cover , G2. 

30 per cent of built-up area with vegetation is added to G2* 

( 3 ) Case 3 : Vegetation cover ,03.^ 

G3 is corrected by subtracting 50 per cent of barren land from G2 > because 
actual barren land does not always involve full of vegetation. 


Table 1 Definition of Vegetation Cover 


^^^Ve^tation Cover 
Class ■ — ■ — 

Gi 

G2 

^3 

( 1 ) Evergreen Forest (Replanted) 

1 

1 

1 

(2) Evergreen Forest (Natural) 

1 

1 

1' 

( 3 ) Deciduous T*orest 

1 

1 

1 

. ( 4 ) Mixed Forest 

1 

1 

1 

( 5 ) Bamboo Forest 

1 

1 

1 

(6) Park 

1 

1 

1 

.( 7 ) Rice field 

1 

1 1 

1 

(8) Crop land 

1 

1 

1 

( 9 ) Orchard 

1 

i 

1 

• ( 10 ) Built-up area with green 

0 

0.3 

0.3 

(11) Built-up area without green 

0 

0 

0 

(12) Street tree 

0 

0 

0 

( 13 ) Water 

0 

0 

0 

( 14 ) Barren land 

1 

1 

0.5 

( 15 ) Disturbed Area 

0 

0 

0 

( 16 ) Other 

0 

0 

0 



MATHER^TICAL MODEL 


The following two models were tested for estimating vegetation cover 
from LANDS AT MSS data, 

(1) Multi-regression model 

Average LANDSAT data with four bands, X (xi , X 2 , X 3 , xi|) in a unit 500 
meters mesh data and the corresponding vegetation cover, y are correlated 
by the multi-regression model as Indicated as follows. 

y = f (X) == aixi+ 32X2+ 33 X 3 + ai|X4+ ao 

Unknown parameters, ao to ai* are estimated by the multi-regression 
analysis, 

(2) Parametric Model 

When variable X, which is derived from LANDS AT MSS data amounts m in 
the unit 500 meters mesh area, let n of m be subject to the condition that 
variable X exceeds parameter k. Supposed that vegetation cover can be 
estimated as n/m, parameter k should be determined so as to minimize the 
standard deviation of residuals between observed vegetation cover G and 
estimated vegetation cover n/m. 

Eight variables were examined in the parametric model, 

Xi = X3 / xi 
X 2 = / Xi 

X 3 - X 3 / X 2 

X^ == X«* / X2 

X5 *= ( X3 + ) / ( Xi + X 2 ) 

Xe === ( XI + X 2 - X 3 - xi^)/( xi + X 2 + X 3 + X 4 ) 

X 7 = aixi + a 2 Xz + asxs + a»fX4 + ao 

( determined by multi-regression model ) 

Xs = hlX* + b2X2 + b3X* + bi*X4 

( second principals determined by principal component analysis) 
(x-x)/a 


' \^ere 



ESTIMATION OF VEGETATION COVER BY MULTI REGRESSION ANALYSIS 


'Multi- regression analysis was applied to 3080 500 meters mesh data with 
three cases of ground truth data. Table 2 summarizes the result of the 
analysis including coefficients ao to and accuracy of estimation represented 
by multi correlation coefficient and standard deviation. Coefficients in the 
lower column represent those for normali 2 ed data. 

Best result was obtained for the case 3, in which multi-correlation 
coefficient y is 0.871 and the standard deviation O is +15.53 per cent. 


Table 2 Result of Multi-regression Analysis 


1 

Coefficient j 

Accuracy 


3-2 

33 

34 

3o 

Y 

a(+%) 

Case 

1 : 

Gi 

-3.844 

-4.839 

3.697 

6.002 

81.368 

0.845 

17.82 




-0.242 

-0.378 

0.331 

0.377 




Case 

2 : 

G 2 

-3.722 

-5.140 

3.874 

5.781 

85.638 

0.856 

16.92 




- 0.237 

-0.406 

0.351 

0.367 




Case 

3 : 

G 3 

-3.340 

-5.229 

3.241 

6.338 

83.417 

0.871 

15.53 




-0.220 

-0.428 

0.304 

0.417 





Fig, 2 shows .comparison between observed vegetation cover and estimated 
vegetation cover for the cross section lines no. 10 , and no. 25 . I'ig- 3 
shows the computer generated map for observed vegetation cover and estimated 
vegetation cover. 


ESTIMATION OF VEGETATION COYER BY FARAMETRIC MODEL 

Most probable parameter k, which would be considered as an criterion to 
decide whether vegetation or not, was determined for each of eight cases of 
variables as mentioned before, so as to minimize the standard deviation 
between observation and estimation. 

Table 3 smumarizes the result of estimation by parametric model. Best 
estimation was obtained from the variable X 3 = X 3 /xa that is the ratio of 
MSS 6 to MSS 5. The standard deviation for this case was +15.36% for the 
most probable parameter k=1.14. This result is of great importance from the 
view point of radiance, that is, since radiances for full scale of CCT data 
in MSS 5 and MSS 6 band correspond to 1.76 and 2.00 mW/cm^.sr respectively, 
the criterion X 3 = X 3 /x 2 jk=l.l 4 coincides with the criterion chat the ratio 
of radiances between the two bands, 1.76 X 3 / 2 . 00 x 2 ^l 5 or X 3 /X 2 S 2 . 00/1 . 76= 
1.14. Therefore, if a slope for the ratio of MSS 6 to KSS 5 with respect to 
radiance is positive, it could be concluded that the corresponding data 
represents vegetation cover. 



Table 3 Kesult of Estimation of Vegetation Cover 
by Parametric Model 



yjase 

l;Gi 

Case 

2 • G 2 j 

Case 

- 3 • G 3 

Variable 

k 

0 

k 

0 

k 

a 

Xi= X 3 /X 1 

0.93 

20.61 

0.90 

21.43 

0.94 

20.05 

X2“ X 4 /X 1 

0.46 

18.97 

0.44 

18.99 

0.45 

18.13 

X 3 = X3/X2 

1.16 

17.27 

1.14 

16.74 

1.14 

15.36 

Xm= Xh/X 2 

0.57 

17.80 

0.55 

17.23 

0.58 

16.12 

Xs= (X 3 +X 4 )/(xi+X 2 > 

0.76 

17.90 

0.75 

17.88 

0.76 

16.66 

Xg" (Xl-hC 2 -X 3 -Xt,) /(X 1 +X 2 +X 3 +X 4 ) 

; 0.14 

17.93 

0.15 

17.86 

0.14 

16.75 

X7= aiXi-ha2X2+a3X3+a4x:ij+ao 

49 

17.58 

47 

17.36 

48 

16.08 

* * * 
Xa- bixi-Hb 2 X 2 +b 3 x 3 + 14 X 4 

- 0.1 

17.85 

; - 0.2 

17.69 

- 0,1 

16.58 


Where X 7 : multi-regression model 

Xe : second principal component 


COMPARISON BETl^EEN MULTI-REGRESSION MODEL AND PARAMETRIC MODEL 

The following comparison between multi-regression model and parametric 
model can be pointed out. 

(1) Multi-regression model utilizes the average data from the stand point of 
statistic, while parametric model accumulates point by point data* 

(2) Parametric model results in a little better estimation, although tKe both 
model have less remarkable difference. Tendencies for error distribution in 
the two models look alike. 

(3) Computing time is much more needed for parametric model, while flexibility 
could be expected more in the parametric model. 


CONCLUSION 

(1) Vegetation cover in Tokyo Districts could be well estimated by both of 
multi-regression model and parametric model. Multi-correlation coefficient 
between observed value and estimated value was 0.87 and standard deviation 
was +^5%. 

(2) Vegetation cover in Tokyo Districts was mapped into five levels with equal 
intervals of 20 per cent. 



(3) Definition of vegetation should be further studied, although the parametric 
model concluded that the decision rule to define vegetation should be whether 
the slope of ratio MSS 6/MSS 5 in the radiance unit is positive or not, 

(4) With detail check for the result of estimation obtained from LANDS AT data, 
actual vegetation cover could be much more represented in the estimated map 
than the ground- surveyed map- 


Line No, 10 



Line "No, 25 


A 



observed value 

estimated value 


Pig. 2 Cross Sections of Vegetation Cover for Comparison 
between Observed Value and Estimated Value 
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FIG. 3 C VEOrrATION CWER ESTimTED BY PAPAIORIC fDEEL 



FIG. 3 B VEGETATION COVER ESTIMATED BY MULTI-REGRESSION ANALYSIS— COLOR PRINTED 
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